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ABSTRACT 

This  document  discusses  methods  for  representing  and  fusing  location 
information  in  combat  models,  distinguishing  between  Monte  Carlo  and  analytic 
models.  The  Kalman  filter  is  emphasized  as  a  practical  and  accurate  method  for 
fusing  new  and  old  information.  One  particularly  useful  Kalman  filter,  the 
Maneuvering  Target  Statistical  Tracker  (MTST),  is  dealt  with  explicitly.  Fused 
information  is  often  represented  graphically  by  ellipses,  the  equivalent  of  a 
Kalman  filter’s  covariance  matrix.  The  connection  between  elliptical 
dimensions  and  kill  probabilities  is  also  reviewed. 

1.  Introduction 

Any  combat  model  needs  to  represent  infonnation  fusion,  whether  implicitly  or 
explicitly.  The  extreme  cases  of  perfect  ignorance  or  perfect  knowledge  are  usually  made 
implicitly.  This  is  the  case,  for  example,  in  systems  of  Lanchester  equations,  or  in 
stochastic  duels.  However,  information  is  sufficiently  important  to  the  Future  Combat 
System  that  combat  models  that  claim  to  represent  FCS  must  deal  explicitly  with 
problems  of  imperfect  information  acquisition  and  fusion.  Hence  the  Army’s  recent 
emphasis  on  realistic  models  of  fusion  within  AWARS  (an  analytic  model)  and 
CASTFOREM  (a  Monte  Carlo  simulation).  This  document  deals  mainly  with  one 
particularly  important  component  of  fusion:  the  fusion  of  information  about  entity 
locations.  The  Monte  Carlo  paradigm  is  particularly  apt  fusion-wise  on  account  of  its 
explicit  representation  of  quantities  such  as  measurements  and  true  entity  locations. 
Fusion  within  analytic  models  is  also  possible,  but  more  difficult. 

2.  Fusion — What  Is  It? 

Fusion  can  be  done  at  multiple  levels  according  to  the  Joint  Director  of  Laboratories 
(Waltz  and  Llinas,  1990).  This  report  is  about  level  1,  the  level  having  to  do  with 
position  and  identity  estimates.  We  suppose  that  here  is  a  collection  of  independent 
“entities”  (squads,  tanks,  command  posts,. . .),  each  of  which  has  a  state  that  is  unknown 
and  may  change  unpredictably  with  time.  Measurements  are  made  on  some  irregular 
schedule,  possibly  in  batches,  and  fusion  is  the  process  of  updating  our  knowledge  of  the 
state  of  each  entity.  It  is  assumed  that  each  measurement  is  caused  by  a  single  entity. 

The  process  of  figuring  out  which  entity  caused  each  measurement  is  data  association. 
Once  the  data  is  associated,  the  process  of  incorporating  it  into  the  state  estimate  for  the 
appropriate  entity  is  filtering.  Both  of  these  real  world  processes  need  to  be  imitated,  to 
some  degree,  in  combat  models.  Although  data  association  must  logically  precede 
filtering,  the  two  parts  will  be  discussed  in  reverse  order. 
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3.  Filtering  of  Position 

An  entity’s  state  vector  may  incorporate  a  diversity  of  properties,  but  the  properties 
nearly  always  include  position  information.  Position  information  can  be  updated  using  a 
Kalman  filter.  Kalman  filtering  is  a  widely  accepted  technique  with  over  four  decades  of 
development  and  experience  behind  it,  and  has  the  further  benefit  of  employing  only 
simple  mathematics,  so  it  is  natural  to  use  it  in  combat  models.  The  basic  idea  in  Kalman 
filtering  is  to  fuse  the  most  recent  measurement  with  the  previous  state  estimate  to 
produce  an  updated  state  estimate.  In  all  cases  the  state  estimate  takes  the  form  of  the 
mean  (p)  and  covariance  matrix  (X)  of  the  distribution  of  the  true  state  of  the  entity.  The 
weight  put  on  the  previous  state  estimate  should  depend  on  how  old  it  is,  as  well  as  on  X, 
since  the  entity  may  have  moved  in  the  meantime.  A  Kalman  filter  thus  needs  two  parts: 
a  motion  model  and  a  measurement  model. 

Motion  model  MTST.  The  Navy  makes  frequent  use  of  a  motion  model  called 
MTST  (Maneuvering  Target  Statistical  Tracker).  MTST  applies  to  targets  moving  freely 
in  a  two-dimensional,  frictional  medium.  The  MTST  state  vector  is  [X,Y,U,V]',  where 
(X,Y)  is  the  entity’s  position  and  (U,V)  is  its  velocity.  Here  the  '  indicates  transpose, 
since  the  state  vector  is  represented  as  a  column  vector. 

Each  velocity  component  is  assumed  to  be  an  Omstein-Uhlenbeck  (OU)  process,  the 
simplest  nonnal,  stationary  process  that  fluctuates  around  0.  An  OU  process  has  two 
parameters:  the  relaxation  time  t  and  the  root-mean-square  velocity  s.  The  average 
velocity  is  0  in  an  OU  process,  befitting  targets  that  move  in  a  frictional  medium  such  as 
land  or  sea,  but  the  root-mean-square  velocity  is  not.  The  relaxation  time  x  is  the  time 
that  it  usually  takes  for  the  velocity  to  change  significantly,  perhaps  an  hour  for  a 
submarine,  0.2  hours  for  a  tank,  or  4  hours  for  a  regiment.  These  two  parameters  (x  for 
time  and  s  for  velocity)  are  the  minimum  number  required  to  provide  scales  for  both  time 
and  velocity.  If  U(t)  is  the  (say  E-W)  velocity  at  time  t,  then  the  OU  motion  model  over  a 
time  perioid  of  length  A  is  U(t  +  A)  =  pU(t)  +  W,  where  p  =  exp(-A/x)  and  IE  is  a  normal 
“process  noise”  with  a  variance  chosen  to  make  the  variance  of  U(t)  be  s~  at  all  times, 
regardless  of  A.  That  variance  is  q  =  s  {  1  -  p")  as  can  be  readily  verified  using  the 
formula  for  the  variance  of  a  sum  of  independent  random  variables.  The  process  noise  W 
represents  the  essential  unpredictability  of  the  entity’s  motion.  The  entity’s  location  X{t) 
is  assumed  to  be  the  time  integral  of  its  velocity.  The  state  update  for  motion  in  a 
Kalman  filter  takes  the  form  p(t  +  A)=Op(t),  and  X(t+A)  =  ®X(t)  ®'  +  Q.  For  MTST, 

®  and  Q  are  the  following  4x4  matrices: 


1  z(\-p ) 

£\  S2 

1  r(l -p) 

and  Q  = 

£\  s2 

P 

s2{\- p2)/ 2 

P 

s2{\- p2)H_ 

The  two  second-order  terms  can  safely  be  taken  to  be  0  when  A  is  small  compared  to  x, 
but  for  completeness,  the  exact  expressions  are 

sx  =  r52(2A-r(3-4/7  +  p2))  and  s2  =  ts2{\-  p)1 . 
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Points  worthy  of  note  about  MTST  are  that: 

•  The  motion  update  algebra  is  simple  to  do  on  a  computer,  amounting  to  one 
exponentiation,  a  few  additions,  and  a  few  multiplications.  The  motion  update  is 
generally  done  only  as  part  of  the  preparation  for  processing  a  measurement,  but 
could  also  be  done  at  the  behest  of  some  person/subroutine/agent  interested  in  the 
current  state  of  the  entity. 

•  When  the  time  since  the  last  measurement  is  small,  x(l-p)  is  approximately  A,  in 
which  case  the  position  updates  for  X  and  Y  take  the  familiar  fonn  of  adding  the 
previous  velocity  estimate  multiplied  by  A. 

•  When  the  time  since  the  last  measurement  is  large,  p  is  approximately  0,  and  the 
position  updates  multiply  previous  velocity  estimates  by  no  more  than  x.  Without 
any  intervening  measurements,  the  location  estimate  moves  in  the  direction  of  the 
most  recent  velocity  estimate,  but  eventually  comes  to  a  stop  as  the  corresponding 
velocity  estimate  approaches  0.  This  is  in  accordance  with  the  OU  process  for 
velocity — regardless  of  the  current  velocity,  the  best  estimate  a  long  time  in  the 
future  (for  an  entity  subject  to  friction)  is  0.  The  phenomenon  has  nonetheless 
caused  some  consternation  among  Navy  officers,  who  initially  expect  the  target  to 
move  indefinitely  at  the  most  recently  estimated  course  and  speed.  The  problem 
disappears  if  there  are  intervening  measurements. 

•  The  two  parameters  cover  stationary  targets  (s=0,  x=  anything),  high  speed  targets 
like  helicopters  (s=50  kt,  x=0.1  hour),  and  medium  speed  entities  like  tanks  (s=5 
kt,  x=.2  hour).  A  separate  logic  is  not  needed  for  every  kind  of  entity.  There 
would  need  to  be  a  database  providing  x  and  s  for  the  various  entity  types.  Since 
it  is  not  essential  to  get  exactly  the  right  movement  model,  it  might  be  possible  to 
use  the  same  default  parameters  for  a  large  class  of  entities  (“tracked”,  say). 

•  The  MTST  code  in  VBA  or  MATLAB  is  available  from  the  author. 

Measurement.  One  of  the  advantages  of  a  Kalman  filter  is  that  there  does  not  have 
to  be  a  single  kind  of  measurement.  The  only  important  feature  is  that  every 
measurement  be  of  some  quantity  or  quantities  corrupted  by  additive,  normally 
distributed  noise.  Possibilities  include: 

•  A  two-dimensional  measurement  of  (X,  Y)  made  by  a  UAV. 

•  A  one-dimensional  measurement  of  arctan  ( Y/X ),  the  bearing  to  the  target,  made 
by  an  ESM  intercept.  This  one  has  been  much  explored  by  the  Navy. 

•  A  one-dimensional  measurement  of  (U2  +  V2)0'5,  the  entity’s  speed,  made  by  a 
seismic  sensor. 

•  A  three-dimensional  measurement  consisting  of  range,  bearing,  and  Doppler  shift, 
made  by  a  radar. 

All  such  measurements  can  be  used  to  update  p  and  X,  thus  fusing  the  measurements  with 
what  was  known  before,  properly  weighted  for  age.  This  all  happens  automatically  in  a 
Kalman  filter,  and  the  mathematics  is  efficient,  once  correctly  implemented.  Details  can 
be  found  in  many  textbooks,  or  in  my  notes  at  http://diana.cs.nps.navy.mil/~washbum/  . 
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4.  Data  Association 

A  measurement  may  contain  clues  or  signatures  that  uniquely  determine  the  target 
that  caused  it.  There  is  no  data  association  problem  in  that  case,  and  fusion  amounts  to 
updating  the  state  of  whatever  entity  was  the  cause.  A  data  association  problem  emerges 
only  when  multiple  entities  might  reasonably  have  caused  the  measurement.  Most  data 
association  schemes  assume  that  the  measurement  to  be  processed  was  caused  by  the 
entity  that  is  in  some  sense  closest  to  it.  In  the  case  of  the  Kalman  filter,  the  best  distance 
to  use  theoretically  is  the  Mahalanobis  distance,  or,  as  I  call  it,  the  dimensionless  shock 
DS.  If  there  are  too  many  possibilities  to  permit  the  calculation  of  DS  for  each  one,  it 
may  be  attractive  to  eliminate  most  of  them  through  some  sort  of  a  gating  scheme. 

Regardless  of  the  association  scheme,  it  is  inevitable  that  an  association  mistake  will 
eventually  be  made.  The  mistake  might  be  to  associate  the  measurement  with  the  wrong 
target,  or  to  fail  to  initiate  a  new  target  track  when  the  cause  of  the  measurement  is  in  fact 
a  new  target,  or  to  fail  to  simply  ignore  the  measurement  when  the  cause  is  not  related  to 
the  military  problem  (a  false  alarm).  Such  mistakes  can  have  very  serious  consequences 
in  the  real  world — one  aircraft  accidentally  shot  down  in  OIF  was  the  victim  of  an 
association  mistake.  The  problem  is  to  imitate  this  in  a  combat  model.  The  method  of 
doing  so  must  depend  on  whether  the  model  is  a  Monte  Carlo  simulation  like 
CASTFOREM,  or  an  analytic  model  like  AWARS. 

5.  Monte  Carlo  Fusion 

One  of  the  beauties  of  dealing  with  a  Monte  Carlo  simulation  is  that  the  true  entity 
causing  any  particular  measurement  is  known  to  the  analyst.  Therefore  one  association 
option  is  to  adopt  the  policy  of  always  associating  data  correctly,  the  “Truth”  policy. 
There  is  much  to  be  said  for  it.  Any  model  of  combat  is  going  to  miss  real  world  aspects 
that  are  important  association  cues,  and  will  therefore  do  a  worse  job  of  association  than 
happens  in  the  real  world.  However,  the  real  world  can  certainly  do  no  better  than  Truth, 
so  Truth  constitutes  an  upper  bound  on  how  well  any  association  scheme  can  do.  There 
are  additional  arguments  that  favor  retaining  Truth  as  one  of  the  association  options.  The 
creation  and  destruction  of  false  target  tracks  would  not  be  necessary,  for  example. 

Using  Truth  does  not  imply  that  the  perception  of  an  entity’s  state  will  always  be 
perfect.  That  perception  depends  on  the  frequency  and  quality  of  measurements,  which 
are  computed  as  usual.  For  entities  tracked  by  a  Kalman  filter,  it  is  worthwhile  to 
compute  DS  for  each  measurement,  as  is  currently  done  in  CASTFOREM.  That  distance 
is  worth  keeping  track  of  even  when  Truth  is  employed.  Statistics  about  DS  will  always 
be  interesting  in  the  sense  that  large  values  correspond  to  mistakes  that  might  have  been 
made,  except  for  the  employment  of  Truth. 

In  addition  to  Truth,  there  should  perhaps  be  some  kind  of  an  automatic  association 
scheme.  The  simplest  such  scheme  would  irrevocably  associate  each  measurement  with 
some  entity,  with  no  possibility  of  subsequent  revision.  Call  this  the  Single  Hypothesis 
Tracker,  or  SHT.  When  the  minimized  DS  is  sufficiently  small  (smaller  than  10,  say), 
the  measurement  would  be  associated  with  the  corresponding  entity.  If  it  is  larger  than 
10,  it  might  simply  be  ignored  on  the  grounds  that,  in  the  real  world,  it  would  probably  be 
a  false  alarm.  In  a  scenario  where  entities  might  plausibly  be  created  or  discovered  as 
time  goes  by,  a  large  DS  might  also  cause  the  creation  of  a  new  entity.  The  SHT  will 
occasionally  make  mistakes,  since  DS  is  occasionally  small  for  a  false  target,  or  large  for 
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a  true  one. 

The  SHT  is  surely  a  simpler  and  less  effective  approach  to  data  association  than  the 
one  proposed  by  Orincon  for  inclusion  in  the  FCS.  As  such,  it  constitutes  a  lower  bound 
on  the  effectiveness  of  the  FCS  system,  just  as  Truth  constitutes  an  upper  bound.  There 
is  thus  a  good  argument  for  including  both  Truth  and  SHT,  especially  if  the  actual  FCS 
system  is  too  complicated  to  retain  in  a  simulation. 

There  are  also  trackers  that  remember  multiple  hypotheses  about  which  data 
association  is  correct,  sometimes  (as  in  Stone  et  ah,  1999)  with  explicit  probabilities  for 
the  correctness  of  the  hypotheses.  I  do  not  see  the  argument  for  including  them  in  an 
existing  combat  model  unless  the  purpose  of  running  the  model  is  to  be  examination  of 
the  effectiveness  of  various  fusion  schemes. 

It  should  be  recognized  that  the  SHT,  or  any  multiple  hypothesis  extension,  will  be  of 
no  use  for  measurements  that  do  not  involve  position.  Kalman  filters  work  only  for  real¬ 
valued  state  components  that  are  plausibly  normally  distributed. 

6.  Analytic  Fusion 

An  analytic  model  such  as  AWARS  (Glasgow  and  Bailey,  2004)  will  have  more 
difficulty  including  association  phenomena  than  will  a  Monte  Carlo  simulation.  The  true 
values  of  state  variables  are  easily  handled  in  Monte  Carlo  simulations  by  sampling  the 
appropriate  distributions.  In  an  analytic  model,  the  only  really  correct  equivalent  is  a 
numeric  integral  over  all  possible  values  of  the  state  variables,  a  much  more  difficult 
operation.  Replacing  distributions  by  their  means  (expected-value  analysis,  or  EVA) 
may  be  the  only  practical  alternative.  In  fact,  EVATruth  is  the  modal  method  for 
representing  fusion  in  analytic  models;  the  association  is  always  correct,  and  actual 
measurements  are  not  represented.  It  is  thus  almost  inevitable  that  an  analytic  model  will 
provide  optimistic  answers  to  problems  where  fusion  is  a  significant  part  of  the  problem. 
The  developers  of  AWARS  recognize  this  (Glasgow  and  Bailey,  2004,  page  2  and 
footnote  5).  It  is  not  practical  to  implement  something  like  the  SHT  in  an  analytic  model, 
so  there  is  no  lower  bound  counterbalance  EVA  Truth. 

There  are  nonetheless  some  powerful  reasons  for  employing  analytic  models,  so 
fusion  needs  to  be  represented  in  spite  of  the  difficulties. 

Kalman  filtering  can  still  be  employed  in  analytic  models,  even  without  specific 
measurements  being  available.  The  best  opportunity  is  with  location  infonnation,  as  in 
Monte  Carlo  models,  and  MTST  remains  an  attractive  candidate  for  the  same  reasons.  A 
Kalman  filter  will  automatically  deal  with  issues  of  data  aging,  can  handle  measurements 
of  multiple  kinds  (bearings,  in  particular),  and  will  automatically  provide  the  kind  of 
information  that  is  required  to  compute  kill  probabilities.  The  state  estimates  themselves 
are  in  some  cases  trivial  because  of  the  EVA  assumption,  but  the  covariance  matrix 
computations  are  still  realistic  and  useful.  In  fact,  it  is  a  curious  characteristic  of  Kalman 
filters  that  the  covariance  computations  are  independent  of  the  actual  measurements,  and 
therefore  of  the  EVA  assumption. 

In  two  dimensions,  the  covariance  matrix  detennines  the  size  and  orientation  of  the 
ellipse  that  graphically  represents  the  target’s  location.  Since  computations  of  kill 
probabilities  rely  on  knowing  the  dimensions  of  that  ellipse,  there  is  a  good  argument  for 
computing  its  dimensions  as  accurately  as  possible. 
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7.  Kill  Probabilities 

There  is  no  point  in  fusing  information  unless  the  information  will  be  used  to  support 
some  decision.  Decision  making  in  the  real  world  will  be  complex,  and  will  involve 
targeting  decisions  about  what  should  be  shot  when  at  whom.  Fusion  levels  far  in  excess 
of  level  1  will  be  involved,  but  level  1  is  still  important  because  it  influences  the  effect  of 
fire.  This  section  deals  with  heuristics  for  allocating  fire  and  assessing  its  effectiveness, 
based  on  the  accuracy  of  target  location. 

Notation.  The  models  below  differ  in  how  lethality  is  modeled,  but  all  share  the 
same  error  structure.  Target  location  error  and  weapon  dispersion  error  (s)  are  each 
assumed  to  be  unbiased  and  elliptically  normally  distributed,  and  all  errors  are  assumed 
to  be  independent.  Target  location  errors  are  assumed  to  be  fused  errors  from  the 
tracking  model,  evaluated  at  the  time  of  fire.  Let 

Gvshorizontal  dispersion  standard  deviation; 
a  ^vertical  dispersion  standard  deviation; 
a  [^horizontal  target  location  standard  deviation;  and 
a  (^vertical  target  location  standard  deviation. 

Note  that  ( U. ,  V)  is  now  the  target  location  error,  rather  than  an  entity’s  velocity.  Two 
separate  two-dimensional  errors  are  involved  because  the  target  location  errors  are 
assumed  to  be  common  to  all  shots  at  a  target,  whereas  the  dispersion  errors  of  different 
shots  at  the  same  target  are  all  independent.  The  assumption  is  that  all  shots  are  taken 
“simultaneously”  in  the  sense  that  no  location  infonnation  is  fused  in  the  time  between 
shots. 

Carleton  weapons.  In  this  case,  Pk=Qxp(-R"/LR~),  where  R  is  the  Euclidean  distance 
between  the  target  and  the  impact  point  of  the  weapon,  and  LR  is  data.  This  formula  is 
also  known  as  the  “diffuse  Gaussian”  formula.  It  is  to  be  distinguished  from  the  cookie- 
cutter  or  “definite  range”  fonnula  where  damage  occurs  if  and  only  if  R  is  smaller  than 
LR.  The  Carleton  expression  generally  results  in  significantly  simpler  expressions  than 
the  cookie-cutter  formula,  especially  when  multiple  shots  are  made  at  the  target 
(Washburn,  2003).  For  a  single  shot  aimed  at  the  center  of  the  target  distribution, 

E(Pk)=LR2/{(LR2+2(ax2+au2))m(LR2+2(aY2+C!V2))V2}. 

The  expected  value  operator  E()  would  normally  not  be  written.  This  is  formula  2-9  from 
(Washburn,  2004),  except  that  LR  =2 b  and  pj*=p  >=0  in  that  document;  i.e.,  the  weapon 
is  aimed  directly  at  the  target.  See  also  chapter  2  of  Eckler  and  Burr  (1972).  In  a  Monte 
Carlo  model  with  imperfect  correlation,  it  is  conceivable  that  one  might  not  be  aiming  at 
the  center  of  the  target  distribution,  either  because  of  a  correlation  mistake  or  because  the 
target  in  question  is  not  the  one  aimed  at.  Whether  the  model  is  Monte  Carlo  or  analytic, 
the  general  form  of  formula  2-9  still  applies  and  is  simple. 

Cookie-cutter  weapons.  In  this  case  there  is  a  simple  formula  only  when  the  total 

horizontal  and  vertical  errors  have  the  same  standard  deviation,  say 

2  2  2  2  2 

g“  =  ax  +cj u  =  or  +CW  .  The  formula  is 

T  R2 

Pk  =  \-  exp( - -) . 

2cr 
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LR  might  be  the  lethal  radius  of  a  weapon,  or  it  might  be  the  radius  of  a  circular 
pattern  of  weapons,  in  which  case  Pk  is  merely  the  probability  that  the  target  is  included 
within  the  pattern.  For  the  ICM  weapon,  this  would  need  to  be  multiplied  by  the 
probability  that  the  target  is  killed,  given  that  it  is  within  the  pattern. 

If  the  horizontal  and  vertical  errors  are  not  equal,  there  is  no  simple  expression  for  Pk 
for  cookie-cutter  weapons. 

Patterns.  When  target  location  errors  are  large,  a  single  round  will  not  be  effective. 
In  that  case  it  may  be  wise  to  use  a  pattern  of  several  rounds,  but  the  subject  is 
complicated  by  the  lack  of  independence  between  rounds  (recall  that  they  all  share  the 
same  target  location  error).  The  weapons  should  generally  be  aimed  in  a  pattern,  which 
raises  the  subject  of  pattern  optimization.  No  combat  model  wants  to  get  into  this  subject 
in  detail,  but  a  simple  formula  relating  the  inputs  to  the  optimized  kill  probability  may 
still  be  useful.  I  suggest  considering  the  “square -root  formula”,  which  is  formula  (3-10) 
in  Washburn  (2002).  The  only  important  quantity  is 

z  =  na/(ln<7u(7V), 

where  n  is  the  number  of  rounds/bomblets  and  a  is  the  lethal  area  per  round.  The 
dispersion  errors  are  unimportant.  The  formula  is 

Pk  =  l-(l  +  >/2z)exp(--s/2z)  • 

If  the  pattern  can  only  be  optimized  to  the  extent  that  all  bomblets  will  be  uniformly 
distributed  over  a  circle  of  the  shooter’s  choice,  then  a  slightly  smaller  kill  probability 
results,  namely 

Pk  =  ( 1  -  exp(- Vz  ))2 . 

Either  of  these  fonnulas  is  simple  enough  to  be  evaluated  many  times  in  the  process 
of  deciding  how  to  shoot  at  a  target,  possibly  by  seeking  the  method  that  is  most  cost- 
effective  in  the  sense  of  maximizing  Pk  per  dollar.  Either  is  sufficiently  accurate  if  the 
number  of  bomblets  is  large  (say,  more  than  10). 

If  exact  pattern  answers  are  required  in  the  Carleton  case  for  10  or  fewer  bomblets, 
see  Washburn  (2003).  The  accompanying  spreadsheet  includes  VBA  code  that  will 
evaluate  Pk  for  an  optimal  pattern,  as  long  as  the  number  of  weapons  involved  is  not 
large. 
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